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Abstract 

The phase of the rotating order parameter in rotating antiferromagnetism is calculated using a combination of mean- 
field theory and Heisenberg equation. This phase shows a linear time dependence, which allows us to interpret 
rotating antiferromagnetism as a synchronized Larmor-like precession of all the spins in the system or as an unusual 
q = (n, 7r) spin-wave around a zero local magnetization. We discuss implications for the pseudogap state of high-Tc 
superconducting materials. Rotating antiferromagnetism has been proposed to model the pseudogap state in these 
materials. 
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1. Introduction 

According to several researchers, the puzzling pseudogap (PG) phenomenon in high-Tc superconductors (HTSC) 
is caused by some sort of hidden order. This is supported by the observation of a depression in the density of states 
at the Fermi level, with no order parameter responsible for this depression observed yet id Hi ■ Rotating anti- 
ferromagnetism (RAF) has been recently proposed as a possible candidate for this hidden order, and several physical 
quantities have already been calculated within the RAF theory (RAFT) with good agreement with available experi- 
mental data [5J, [ZD . RAF is one of several other proposals for the PG (see Ref. Q4£| for a discussion). Contrary to 



theories of circulating currents 081 19[ [lOfl, RAF is based on the concept of an order parameter that has a finite magni- 
tude below a critical temperature but a time-dependent phase ll ill . Note that all the physical quantities that have so far 
been calculated within RAFT do not depend on the phase of the order parameter in RAF [5,6, [H E Ellil. The 
lack of the time dependence profile for this phase limited however the full understanding of the nature of RAF. The 
purpose of this work is to calculate this phase as a function of time using a combination of RAFT and the Heisenberg 
equation. We show that it varies linearly with time. As a consequence of this time dependence, RAF can be interpreted 
as a (n, n) unusual spin wave around a zero local magnetization or as a synchronized Larmor-like precession of all the 
spins in the system. Because the phase of this order parameter is time dependent, it was not possible to calculate it in 
RAFT alone, which is a mean-field approach. 

This paper is organized as follows. First in Sec. 12. II we rederive RAFT using the spin ladder operators, which 
are necessary for the phase calculation. In Sec. 12.21 we review RAFT. Then in Sec. 12.31 we use the Heisenberg 
equation to get the time dependence for the spin ladder operators, which yields the time dependence of the phase 
of the rotating order parameter. In Sec. 12.41 the interpretation of RAF as an unusual (n, n) spin wave is explained. 
Finally, conclusions are drawn in Sec. [3] 

2. Approach 

As we are only interested in understanding the nature of the PG phase of HTSCs in this work, we restrict ourselves 
to the non superconducting phase. Consider the t-t' Hubbard model in two dimensions: 

H = - t Yj c l°- c i> cr ~ t ' Yj • |KC - 

(i,j)<r «/J»£r 
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where (i, j) and ((/, j)) designate summation over nearest and second-nearest neighboring sites, respectively, t and 
t' are electron hopping energies to nearest and second-nearest neighbors, respectively. Because the phase of RAF is 
related to the spin ladder operators, it is useful to rewrite Hamiltonian ([TJ using these operators. 

2.1. Rewriting the Hamiltonian using the spin ladder operators 

Using the spin ladder operator written in second quantization St = cT^c,-^, the onsite Coulomb repulsion term 
Un^nii canon one hand be cast in the form Un^nn - Un^— US tSj and on the other hand as Un^n^ - Un^—US'St . 
Summing and dividing by 2 yields the symmetrized expression Un^n^ - ^(h/t + n,x) - j(StSj + SjSt). The latter 
can be proved by calculating the action of each side of the equality on the possible states {|0), | T)> | |), | TDK an d 
noting that StSj\ fl) = StS~!~\ tl> = due to Pauli exclusion principle, and S^IO) = c^c^O) = 0. For our many- 
body system, sites are neither full nor empty, but are on average occupied by a density smaller than 1 away from half 
filling. Therefore, the terms S^St and S^St, which are responsible for onsite spin-flip excitations, will contribute 
by lowering energy for the sites that are partially occupied by the same density of spin up and down electrons. One 
can decouple this term in mean-field theory using (Sj) = (cT.C;^), which leads to a collective behavior for the spin- 
flips, and the results obtained in this way are the same as in RAFT IdI la 171 Il2l Il3l 1 1-411 . In this state, a spin flip 
process at site i is simultaneously accompanied by another one at another site j; the occurrence of the spin flips is 
synchronized. Thermal motion has obviously an effect on this order as it does on conventional orders; i.e., above a 
critical temperature (identified with the PG temperature) the spin-flip processes become uncorrelated, leading to the 
disappearance of the long-range non conventional order. The spin-flip processes, which are purely quantum, continue 
to exist even above this critical temperature, but in an incoherent disordered manner. The occurrence in RAFT of a 
second-order phase transition at the PG temperature is consistent with experimental data supporting its existence fl. 

2.2. Review of RAFT 

We rederive RAFT, which deals with the static part (magnitude) of the order parameter (Sf), using the spin ladder 
operators then for the dynamic (phase) part we will use the Heisenberg equation to find its time dependence. To 
the best of our knowledge the combination of mean-field theory and the Heisenberg equation of quantum mechanics 
constitutes a novel approach for the PG in HTSCs. 

The parameter Q; = (c^fcT,) = -(St) = |<2|e i( *' is defined in order to carry on a mean-field decoupling of the t-f 
Hubbard model. Consider the ansatz where 0,- - (pj = n, with i and j labeling any two adjacent lattice sites. Except 
for this difference of n between the phases of the order parameter on two adjacent sites, the phases 0, = <p are site 
independent and assume any value in [0, 2n]. The normal state Hamiltonian in RAFT JUSO] is 

ff« V ^"WPt + NU Q 2 - NUn 2 , (2) 

keRBZ 

where is the number of sites, and n = <n,>) is the expectation value of the number operator. Because of an- 
tiferromagnetic correlations the lattice consists of two sublattices A and B, even though there is no long-range 
static antiferromagnetic order. The summation runs over the reduced Brillouin zone (RBZ). The Nambu spinor is 
*P k = (c^ c^-j c^j ), and the Hamiltonian matrix is 

e Qe^ ^ 
-H' -Qe^ 
-n' e 
Qe-* e -y! ) 

yielding the energy spectra E ± (k) - - / u'(k)±£' 9 (k), where yu'(k) = fi-Un+4t' cos k x cos k y , E q (k) = y/e 2 (k) + (UQ) 2 , 
and e(k) = -2f(cos k x + cos k y ). Because the energy spectra E ± (k) do not depend on the phase <p one should be able to 
transform H to a matrix that does not depend on the phase. This can indeed be done using the spin-dependent gauge 
transformation c, j — > e"^ 2 c,- ^ and c, j — > e~'^ 2 c, j. This transformation is equivalent to performing a rotation by angle 
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-tf> about the z axis for the x and y components of the spin operator. Indeed, upon using this gauge transformation, the 
spin ladder operators transform according to St — > e _!<s St and S7 — > e^Sj, which yields: 

S* \ / cos</> sin0 \/ S? 
S 31 I y -sin</> cos0 l\ S? 

The thermal averages of Sf and S ■ are given by 

(Sf) (S)) 

= Qcoscf), = -Qsin<f>, i e A, or 

(Sf) (S y ) 

— l — = -Qcos<f), — '- = Qsincf), ieB. (3) 

h h 

Note that (S z ) = for i in both sublattices. Because the phase <p assumes any value between and 2n (see below), 
rotational symmetry will not look broken for times greater than the period of rotation. However if the typical time 
scale of a probe is much smaller than this period symmetry may appear broken. 

2.3. Calculation of the time dependence of phase <p 

The magnitude Q, which was calculated using the minimization of the mean-field free energy JUSO], behaves as 
in a second-order phase transition in agreement with experimental evidence in fl. 

Next we calculate the phase using the Heisenberg equation = ^[St,if]. We consider the limit where electron 
hopping is neglected in comparison to ~ (S + S 7 +S tS + ). The limit considered here is U ~ 3?-5f; this is an intermediate 

^ J J J J 

coupling limit where U > t but smaller than the bandwidth ~ 8f when r <K t. It is justified to use this approximation 
because spin dynamics is faster than charge dynamics; i.e., an onsite spin flip needs a time r ~ h/U to be realized, 
while a charge hopping between adjacent sites takes a longer time t ~ H/t, (U > t). In the Heisenberg equation 
the undecoupled interaction is used instead of RAFT's Hamiltonian (0 in order to treat as best as possible quantum 
fluctuations. To carry on the calculation, we keep in mind that any site j is on average only partially occupied, and 
that \{SJ)\ < h/2. For this reason, terms like S^S^Sj and S^SjS^ should be kept until the end (these terms normally 
give zero when acting on a spin up state, but a nonzero contribution is expected when applied to a partially occupied 
state where thermal averages are meaningful and suitable). In the commutator of the Heisenberg equation [S^,H] * 
-f [St, (StSj+SjSt)], we need to calculate [St, (S+St +SjSt)] = [St,S+Sj] + [S+,SjSt] = 2h{S ]S\-S : ,S ]). 
Using the fundamental commutation relation [S*, S t] = HS t, one gets S tS z . + S Z jS t = HS t + 2S tS z -, which leads to 

dS+ Ijj 2U \ 



dr 

Again we stress that this equation is obtained in the intermediate coupling limit (U smaller than the bandwidth but 
higher than hopping energies), where spin dynamics is not governed by the Heisenberg exchange coupling ~ t 2 /U 
suitable for the strong coupling limit. Eq. © gives zero when acting on state | f ) or ] |). However, for a collective 
state where any site is only partially occupied, one has to take the thermal average of Eq. (|4). One then replaces S z 
by its RAFT's thermal average, which is zero. Integrating Eq. gives for the thermal average 

(St( T )> * (S+(0))e iUT/H , (5) 

which yields <p = Ur/hmodulo 2n when (S t(0)> is identified with |<S t( T ))|, (-RS *(r)>|), for sublattice A, (B), and e^ 
with e' UT l h . The angular frequency is thus 0) s f = U/h, and period T s j = 2jtH/U is the time required to perform a spin- 
flip process, or the time needed for the rotating order parameter (S f-^) to complete a 2n revolution in a classical point 
of view. The magnetic configuration Q takes on the following form (S x )/H = Qcos(a> s fT), (S y -)/H = -Qsin(o) s fT) 
for i in sublattice A or (S x )/h = -Qcos(u) s fT), (S y -)/H = Qsm{oj s fT) for i in sublattice B, and (S?) = for i in 
sublattice A or B. 
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2.4. Interpretation of rotating antiferromagnetism 

To illustrate well rotating magnetic order, consider first the much simpler example of the time evolution of a single 
spin in a magnetic field B along the z-axis, with the initial state given for a spin pointing in the positive ^-direction 
by \S x , +) = + I - ))- The time-dependent expectation values of the spin components are (S x ) = jcos(tot), 

(S y ) = § sin(wf), and (S z ) = 0, with to = ^. e and m e are the charge and mass of the electron, respectively, and c 
is the speed of light. Classically speaking, the spin is confined to rotate about the z-axis in the xy plane with Larmor 
angular frequency a>. A rotating ferromagnetic state can be realized by placing N such states with the same frequency 
on a lattice made of N sites. For a rotating antiferromagetic state, opposite initial states (±|5^, +): spins point in 
opposite directions on the x-axis) are required on each two adjacent sites of the lattice. To relate RAF to spin flip 
processes, we note that (S ± ) = (S x ) + i{S y ) = je ±uo ' in this example. Note that in this example is model independent, 
which may indicate that all model parameters will do is changing multiplying physical factors, not the physics itself. 
In a given model, a coupling is necessary for providing the building bloc for RAF, which is the precession of a spin 
(with no local magnetization) for each lattice site. The RAF state constructed in this way shows a hidden order that 
can be realized even at finite temperature without violating the Mermin-Wagner theorem [16]. 

The example above allows us to interpret RAF as a state where spins precess collectively in a synchronized 
manner in the spins' xy plane around an effective staggered magnetic field B = m e cU/h\e\ caused by onsite Coulomb 
repulsion. For our many -body system, h/2 in (S ± ) = |e ±,a " is replaced by the magnitude of the RAF order parameter 
Q, which can assume values smaller than 1/2 due to thermal averaging. This state is strongly doping dependent. When 
doping increases, Q rapidly decreases then vanishes at a doping identified as the quantum critical point underneath 



the superconducting dome |5|, |a |7|, [12J, |13j, [14j . In comparison to ordinary spin waves in an antiferromagnet, RAF's 
state could be viewed as a q = (n, n) spin wave in an antiferromagnet with zero magnetization. Note however that for 
our system (where (Sp = 0), spin-wave theory is not applicable because the spin-wave theory is built around a stable 
nonzero (Sp state. 



3. conclusion 

The rotating antiferromagnetism theory and Heisenberg equation are combined in order to calculate the phase of 
the rotating order parameter. This phase behaves linearly in time. This allows us to interpret rotating antiferromag- 
netism in terms of a Larmor-like spin precession about an effective magnetic field, which is proportional to onsite 
Coulomb repulsion. Another way to see rotating antiferromagnetism is as an unusual spin-wave at q = (n, n) around a 
zero magnetization. This work was necessary for unveiling the nature of rotating antiferromagnetism, which has been 
proposed for explaining the pseudogap behavior in high- Tc materials. Rotating antiferromagnetic order is an example 
of hidden order, which is a serious candidate for the PG state in HTSCs. This is supported by the good success of 
the rotating antiferromagnetism theory in the calculation of thermodynamics lHH^Il optical conductivity iQj, 151. 
Raman fl3fl . and angle-resolved photoemission spectroscopy properties J3- 

Author wishes to thank A.-M. S. Tremblay for helpful comments on the manuscript. 
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